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We study the octet baryon electromagnetic properties by applying the covariant spectator quark 
model, and provide covariant parametrization that can be used to study baryon electromagnetic 
reactions. While we use the lattice QCD data in the large pion mass regime (small pion cloud 
effects) to determine the parameters of the model in the valence quark sector, we use the nucleon 
physical and octet baryon magnetic moment data to parametrize the pion cloud contributions. The 
valence quark contributions for the octet baryon electromagnetic form factors are estimated by 
extrapolating the lattice parametrization in the large pion mass regime to the physical regime. As 
for the pion cloud contributions, we parametrize them in a covariant, phenomenological manner, 
combined with SU(3) symmetry. We also discuss the impact of the pion cloud effects on the octet 
baryon electromagnetic form factors and their radii. 
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I. INTRODUCTION 

Low-lying baryons are classified into the octet 
(spin 1/2) and decuplet (spin 3/2) baryon members. 
This is based on SU(3) symmetry, which is presently un- 
derstood in terms of quantum chromodynamincs (QCD). 
To study the electromagnetic structure of the light 
baryons, they are probed by electron beams, and the form 
factors are measured as functions of momentum trans- 
fer squared q^ = — Q^, and they reveal the nonpoint- 
like structure of the baryons. In the octet and decuplet 
baryons, only the weak and electromagnetic structure of 
the nucleon has been well measured in experiments at 
finite Q^ . Our present knowledge of the electromagnetic 
structure for the octet and decuplet baryons is restricted 
to magnetic moments of some octet [H, l2l| and decuplet (3| 
baryons. Thus, except for the nucleon system [1|, studies 
of the weak and electromagnetic structure for the other 
octet baryon members are very scarce |4|-[l6|. Therefore, 
more data and studies of the octet baryon structure are 
desired. 

In this work we study the electromagnetic structure 
of the octet baryons using the covariant spectator quark 
model [il-y|. In particular, one of our goals is to ob- 
tain covariant parametrization associated with the va- 
lence quark degrees of freedom, that provides insight into 
the octet baryon electromagnetic internal structure. This 
will open tremendous possibilities for future applications. 
A study of the decuplet baryons with the same formalism 
was already performed successfully in Ref . [3[ . 

Generally, phenomenological treatment of baryon 
structure based solely on the valence quark degrees of 
freedom, can be improved by the inclusion of the meson 
cloud effects, guided by chiral effective field theory [l7| . 
The effects of the meson cloud are particularly important 
for the neutron electric form factor fill, and the electro- 



magnetic transition, -fN -^ A(1232) [I8|-l20|. Although 
chiral perturbation theory is very useful to infer the Q^ 
dependence of the nucleon form factors, it is usually ap- 
plicable in the very low Q^ region (< 0.4 GeV^), and can- 
not be used to define the scale of separation between the 
valence quark dominant region from that of the meson- 
baryon excitations (meson cloud effects) [21|- The sepa- 
ration is possible only within a specific model. QCD in 
the limit of high Q^ , and many phenomenological calcula- 
tions suggests that the meson cloud effects should falloff 
with increasing Q^, leading to the valence quark domi- 
nance in the high Q^ region. However, the magnitude, 
sign, and the rate of the falloff of the meson cloud excita- 
tions depend on models [g, l2l| - [33 |. Among the mesons, 
the pion is expected to give the most important contri- 
butions according to chiral perturbation theory [21- 23J. 
For example, the importance of the effects in nucleon 
magnetic moments in different models can be found in 
Ref. 2]. However, for the nucleon, an accurate descrip- 
tion of the form factors can also be achieved by a model 
without the pion cloud in the physical and the lattice 
QCD regimes P,[3l], e.g., model II in Ref. ^. This sug- 
gests that the effects of the pion cloud can possibly be 
small for the nucleon form factors in the low Q^ region. 
Although the pion cloud effects cannot be extracted di- 
rectly from the experimental data, we assume that the 
octet baryon systems can be described by a mixture of 
the valence quarks and the pion cloud, and use lattice 
QCD data to constrain the valence quark contributions, 
and then to extract the pion cloud effects. 

In this study we use the covariant spectator quark 
model which is inspired by the covariant spectator the- 
ory [3J|. The covariant spectator quark model has 
been successfully applied for studying the electrornag- 
netic properties of several baryons [l|-[3, [l8l - l20l . [33l ISSl - 
|44| . A baryon in the model is described by a wave 
function for the quark-diquark system parametrized in 



a covariant manner. The photon-quark couphng is de- 
scribed based on vector meson dominance (VMD) for a 
constituent quark. Because of this feature, the model can 
be extended easily to the lattice QCD regime 0, |20, [33| , 
as will be explained later. Starting by a simplified model 
for the nucleon based on an S-state configuration for the 
quark-diquark system, we extend the model for the octet 
baryons first, and next to the lattice QCD regime. That 
procedure allows us to use the lattice QCD simulation 
data with large pion masses to constrain better the pa- 
rameters of the model, instead of using only the phys- 
ical data. The final parametrization for the wave func- 
tions obtained in the lattice regime, can then be extrapo- 
lated to obtain the physical octet baryon wave functions. 
This is particularly useful, since there are no experimen- 
tal data for E, A and S baryons for Q^ > 0. However, 
this procedure can only provide the contributions from 
the valence quarks, since still lattice QCD simulations are 
presently performed with large pion masses, and thus the 
effects of the pion cloud are suppressed [43| . To take into 
account the effects of the pion cloud which are very im- 
portant for the physical regime, we first use a simple, co- 
variant parametrization for the pion cloud contributions 
for the nucleon electromagnetic form factors. Then, the 
parametrization for the nucleon is extended to the other 
octet baryon members using SU(3) symmetry. 

This article is organized as follows. In Sec. |TT] we de- 
scribe the electromagnetic current, and discuss the sepa- 
ration of the photon couplings with the quarks, and those 
with the pion cloud. In Sec. HII] we give detailed expla- 
nations of the covariant spectator quark model extended 
for the octet baryons, both for the physical and lattice 
QCD regimes. In Sec. IIVI we discuss the pion cloud 
contributions, and give phenomenological parametriza- 
tion for the pion cloud contributions. In Sec. FVl explicit 
parametrization for the octet baryon wave functions and 
pion cloud contributions for the electromagnetic form fac- 
tors are presented. The results for the bare and dressed 
electromagnetic form factors of the octet baryons are pre- 
sented in Sec. lVIl Summary and conclusions are given in 
Sec.lvni 



II. VALENCE QUARK AND PION CLOUD 
CONTRIBUTIONS FOR THE CURRENT 

The current associated with the elastic electromagnetic 
interaction with a baryon B, with spin 1/2 positive parity 
and mass M^, can be represented in general as 



J'B=FlB{Q^)r 



F2b{Q 



a.^o-^^q^ 



2Mf 



(1) 



where Fib and F2B are respectively the Dirac and Pauli 
form factors which are the functions of Q^ — —q^, where 
q = P_|_ — P_ , with P_|_ (P_ ) being the final (initial) mo- 
mentum. Omitted in Eq. ([ij are the initial and final 
state Dirac spinors function of P± and the spin projec- 
tions. For simplicity we represent the current in units 
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FIG. 1: Electromagnetic interaction with the baryon B within the 
one-pion loop level (pion cloud) through the intermediate baryon 
states B' . A diagram including a contact vertex fizBB', as de- 
scribed in Ref. Ql , is not represented explicitly, since the isospin 
structure is the same as diagram (a). See Ref. @| for details. 



e = •\/47ra, with a ~ 1/137, the electromagnetic fine 
structure constant. At Q^ — 0, these form factors are 
defined as 



Fis(0) = eB, F2b(0) = kb, 



(2) 



where bb is the baryon charge in units of e and kb 
is the baryon anomalous magnetic moment in natu- 



ral units 



2Me 



An alternative representation of the 



electromagnetic form factors of the baryon B is the 
Sachs parametrization in terms of the electric charge 
Geb = Fib ~ jjjtF-zb and magnetic dipole Gmb = 
Fib + F2B form factors. 

In a quark model the electromagnetic interaction with 
a baryon B may be decomposed into the photon interac- 
tion with valence quarks, and with sea quarks (polarized 
quark- ant iquark pairs). The latter can be interpreted 
as virtual mesons which dress the baryon valence quark 
core. The photon couplings with the intermediate meson- 
baryon states can be described by effective field theories 
that are perturbative for the low meson energies and mo- 
menta. According to chiral perturbation theory, the most 
important meson for a given reaction is the lightest one, 
the pion. This is also expected to be the case for the octet 
baryons. Thus, we can describe the electromagnetic in- 
teraction for a member B of the octet baryons using a 
current. 






Zb ^b + J^ + J'^b 



(3) 



where J^^ stands for the electromagnetic interaction 
with the quark core without the pion cloud, and the re- 
maining terms are the interaction with the intermediate 
pion-baryon (ttB) states. (See Fig. [T]) In particular, J^ 
represents the direct interaction with the pion, and J^^ 
the interaction with the baryon while one pion is in the 
air. The factor Zb is a renormalization constant, which 
is common to each isomultiplet: nucleon [N), S, A, and 
S. Zb is related with the derivative of the baryon self- 
energy [2]- 

In an additive constituent quark model, the current 
J^g is given by the sum of the individual quark current. 



The electromagnetic interaction processes for the nucleon 
in the covariant spectator quark model 1* is presented in 
Fig. [21 The decomposition given by Eq. ([3]) is justified 
when the pion is created by the overall baryon, but not 
by a single quark. The processes where pious are cre- 
ated and absorbed by the same quark are included in the 
constituent quark internal structure, and thus included 
in the current J^g. 




III. SPECTATOR QUARK MODEL FOR THE 
OCTET BARYONS 

In the covariant spectator quark model a baryon B is 
described as a system with an off-mass-shell quark, free 
to interact with photons, and two on-mass-shell quarks. 
Integrating over the two on-mass-shell quark momenta, 
we represent the quark pair as an on-mass-shell diquark 
with an effective mass tti^i, and the baryon as a quark- 
diquark system [l|. This quark-diquark system is then 
described by a transition vertex between the three-quark 
bound state and the quark-diquark state, that describes 
effectively the confinement [l|, Q . 



A. Octet baryon wave functions 



FIG. 2: Electromagnetic interaction with the baryon B in 
a impulse approximation. P+ (P-) represents the final (ini- 
tial) baryon momentum and k the momentum of the on-shell 
diquark. The baryon wave function is represented by V^s. 
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ub{P,s). (10) 



In the above ub{P,s) is the Dirac spinor of the octet 
baryon B with momentum P, spin s, and Sx{P) the di- 
quark polarization in the fixed-axis representation |ll.l35|. 
In Ref. PI it is shown how Eq. ^ generalizes the non- 
relativistic spin wave functions. 



The simplest representation for a quark-diquark sys- 
tem with spin 1/2 and positive parity is the S-wave con- 
figuration. The wave function for an octet baryon B with 
momentum P and the internal diquark momentum fc, can 
be represented in the S-state approximation fl', '2] , 

^b{P, fc) = ^ {^°s\Ma} + (l>'s\Ms}} ^/jb{P, k), (4) 

where \Ma) and \Ms) are the flavor antisymmetric and 
mixed symmetric states, respectively, and (j)g are the 
spin (0 and 1) wave functions. ^'BiP^k) is a scalar func- 
tion of P and k, and it reflects the momentum distribu- 
tion of the quark-diquark system. 

Explicit baryon flavor wave functions are presented in 
Table HI The spin wave functions are given by 1] 
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4(+5) = -^[(n + ;t)t-2m], 
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This nonrelativistic structure is generalized in the co- 
variant spectator quark model [l| as 



ii 



= ub{P,s), 

= -eT{P)U%{P,s), 



(9) 



B. Electromagnetic current 

Taking into account that the wave function '^ b is writ- 
ten in terms of the wave functions of a quark pair (12) 
and a single quark (3), one can write the electromag- 
netic current associated with the baryon i? in a impulse 
approximation pj, i3| , 



'-'OB 



^ 3J2 f '^B{P+,k)jj^^B{P-,k), (11) 

T^ J k 



where j^ is the quark current operator, P+ {P-) is the 
final (initial) baryon momentum and k the momentum 
of the on-shell diquark, and F = {s, A} labels the scalar 
diquark and the vectorial diquark polarization A = 0, ±. 
The factor 3 in Eq. ([TT|) takes into account the contribu- 
tions for the current from the pairs (13) and (23), where 
each pair has the identical contribution with that of the 
pair (12). The polarization indices are suppressed for 
simplicity. The integral symbol represents 
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where Ed — yr 

Generally, the baryon electromagnetic current pT|) can 
be expressed as 
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where Sob and kqb are the functions of Q^ , and respec- 
tively correspond to the valence quark contributions for 
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TABLE I: Flavor wave functions of the octet baryons. 



the Fib{Q^) and F2b{Q^) form factors. To represent 
these quantities for Q^ = 0, we suppress the tildes. Note 
that in Eq. (J13p we omit the baryon spinors as in Eq. ([Ij . 



C. Quark current 



The quark current operator j^ has a generic structure, 



j'i^nii' 
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2Af, 
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where Mm is the nucleon mass and ji (i = 1, 2) are SU(3) 
flavor operators acting on the third quark of the [Ma] or 
[Ms] state. In the first term 4q^ /q^ is included for com- 
pleteness, but does not contribute for elastic reactions. 

The quark current ji {i = 1,2) in Eq. 1^^, can be 
decomposed as the sum of operators acting on quark 3 
in SU(3) flavor space 



ji — Tifi+^O + ofi-^i + fi/joAs, 



(15) 



where 



An = 




A3 = 




(16) 



are the flavor operators. These operators act on the 
quark wave function in flavor space, q = (uds)"^. 

The functions fi±{Q'^) {i — 1,2) are normalized by 
/i„(0) = 1 (n = 0,±), /2±(0) = K±, and /2o(0) = k,. 
The isoscalar (k+) and isovector (k_) anomalous mag- 
netic moments are defined in terms of the u and d 
quark anomalous magnetic moments, k+ = 2k„ — k^ and 
K- = iKu + -^Kd- In the previous works the quark anoma- 
lous magnetic moments were adjusted to reproduce the 
experimental magnetic moments of the nucleon and the 



ri~ [ij, y]. In this work however, we will readjust the 
u and d quark anomalous magnetic moments as will be 
explained later. 

To see explicitly the quark flavor contributions for the 
electromagnetic current (|14p , we sum over the quark fla- 
vors following Refs. [2, 3], and get the coefficients 






(17) 
(18) 



for i — 1,2. The results, corresponding to the states 
given in Table HI are presented in Table [Til 



D. Valence quark contributions for the 
electromagnetic form factors 

Using the expressions derived in the previous work 
for the nucleon form factors in the S-state approach [l], 
we obtain the corresponding expressions for the octet 
baryons B by replacing the nucleon coefficients jf and 
Ji (i = 1, 2) by the respective baryon state 
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the overlap integral between the initial and final scalar 
wave functions. The normalization of the wave function 
leads to B{0) = 1. The expressions in Eqs. ([HI) and ((20)) 
given for the nucleon pj are briefly reviewed in the Ap- 
pendix. 
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and 77, can be described by the two independent coeffi- 
cients D and F based on SU(3) symmetry [Uliil. The 
couphng constant of the pion (tt) and baryons (B and 
B'), gwBB', can be represented in terms of the ratio, 
a = -p^To ^^"^ ^ global coupling constant g = gy^NN: 
the ttNN coupling constant. As a result, the interaction 
currents and baryon self-energies, at the one-pion loop 
level, which depend on the coupling constants of, ttNN, 
ttAS, ttSS and ttSS, can be expressed in terms of the 
independent coefficients /3b for B = A,T,, S, 



TABLE II: Mixed symmetric and antisymmetric coefficients 
for the octet baryons appearing in Eqs. (|17[) and <\lS\i. 



Another, possible electromagnetic transition between 
the octet baryon members is 7* A — ^ E°. This transi- 
tion is very interesting, and will be studied in a separate 
work [4^, since it is an inelastic reaction. 

The expressions for cqb and kqb [Eqs. (ITOl) and d^ ] 
are particularly simplified for Q^ = 



eoB 
kqb 



[J^+J^) 
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(22) 



We require that the bare charge eos and the dressed 
charge es are the same 



eos = e_B. 



(23) 



To get the numerical results, we must specify the func- 
tions f^n{Q'^) (i = 1,2, n = 0,±) and iJBiP,k). The 
explicit expressions will be given in Sec. |Vl 

The form factors described in this section, correspond- 
ing to Eqs. ([T9|) and ([20]) . include only the valence quark 
contributions. For a realistic estimate we need to include 
the pion cloud effects explicitly. 



IV. PION CLOUD CONTRIBUTIONS FOR THE 
ELECTROMAGNETIC FORM FACTORS 

We discuss here the pion cloud contributions for the 
electromagnetic current and form factors, represented by 
the diagrams in Fig. [TJ Following Ref. p| , we assume the 
pion as the dominant meson excitation to be included in 
the octet baryon form factors. Then, the meson cloud 
contributions for the octet baryon electromagnetic form 
factors can be described in terms of 6 independent func- 
tions of Q^, related to the pion- baryon Feynman integral, 
as will be described next. 



/?s = 4(l-a)2, 
;3h = (1 - 2af. 



(24) 

(25) 
(26) 



The factor for the nucleon is (3n = 1- Absorbing the 
global coupling constant g in the functions associated 
with the pion-baryon loops, we can express the pion cloud 
contributions entirely in terms of a. 

In the exact SU(3) symmetry limit all the octet baryon 
masses are the same. We break the symmetry by the 
pion cloud effects on the mass, but only by the one-pion 
loop in the self-energies. Following Ref. [2|, we represent 
the mass of the octet baryon B as AIb = Mq + Yj{Mb), 
where Mq is a mass parameter and S(Mb) the self-energy 
at the physical baryon mass point. Then, we can write 
S(Mb) = GobBq, with Bq being a scalar integral and 
Gob a factor that includes the couplings of the pion with 
the baryon (see Table 2 in Ref. [4]). Using the value, 
a = 0.6 given by an SU(6) quark model [2J| , we get Mq = 
1.342 GeV and Bq = -0.127 GeV, which can describe the 
octet baryon masses within an accuracy of 7%. 

The choice of a = 0.6 given by SU(6) symmetry, 
defines the strength of the pion cloud contributions as 
/3e = if, /3a = i and /3h = 5^ [1, |2|. Note that, the 
intermediate baryons in the one-pion loop diagrams in 
Fig. [TJ are always those in the same isomultiplet with 
the external baryons (iV, S or S), with the exception of 
the A-E mixture for the E case, where the mass difference 
of them is small. 

A comment on the value of a is in order. Contrary 
to the previous study of the spectator quark model p| , 
where a ~ 0.69 was obtained by fitting to the octet 
baryon masses and octet baryon magnetic moments, here 
we use the value obtained by SU(6) symmetry. This en- 
ables us to treat all the octet and decuplet baryons in a 
unified manner based on SU(6) symmetry. 



A. Pion-baryon couplings in SU(3) 



B. Pion cloud dressing 



The meson-baryon interaction vertices between the 
baryon octet and the pseudoscalar meson octet, tt, K, K 



The pion cloud corrections, namely the coupling of the 
photon to the pion J^, and the coupling to the interme- 



diate baryons J^g, can be written [2| 



J^ = [BlJ^ + B; 
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(27) 



(28) 



In the above, Bi, Ci and Di (i — 1, 2) are arbitrary func- 
tions of Q^ (as Gob and kqb are), and Gt^b^ G^b, G^b 
are the coefhcients that depend on the baryon flavors 
(B = A^, E,A,S). We assume that the functions Bi^ Ci 
and Di are only weakly dependent on the baryon masses, 
and the same for all the octet baryons as in Ref. [2| . That 
allows a description of the pion cloud dressing with a re- 
duced number of coefficients. We write Bi, d, and Di 
to represent, respectively, the functions Bi, Ci and Di 
at Q^ = 0. The coefficients Bi and B2 are proportional 
to the pion electromagnetic form factor F.^{Q'^), but we 
absorb it in the definition of Bi for simplicity. 
Combining Eqs. 0, (I2Z1) and ^, we get 



J'^B = ZB{eoB + G^BBi + GeBCi + G^BDl]l' 



Zb \ Hob + G-^bB2 + GesC'2 + GkbD2 



2Mb ■ 
(29) 



The renormalization constant Zb is determined by the 
relation between the dressed form factor Fib a.t Q^ — 
and the charge e_B [^is(O) = es], or by the self-energy. 
(See Ref. [2j for details.) The condition for the nucleon 
leads to 



Di=0, 



Bi=C, 



(30) 
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TABLE III: Coefficients Crs and GzB, where z = e,K for the 
octet baryons _B = A^, A, S and H. See also Ref. Q] for details. 



convenient to introduce a general operator decomposition 
for the bare current. We use zb to represent cb or kb- 
For the N and S isospin doublets, we use the standard 
isoscalar-isovector notation. 



ZB = UzB+ZB-ra) 



(33) 



At Q^ = we have e|^ = e^ = 1, and e| — — e| — —1. 
The A case is given by the scalar functions cqa and kqAj 
with CA = eoA = 0. 

For the E isospin operators, the decomposition of the 
three states (0, ±) can be given by 



z^ = 41 



I'ziJa + z^J 



2 t2^ 



(34) 



where J3 is the third component of the isospin 1 operator. 
With this notation we get 



and gives 



Zat = 1 + 3Bi 



(31) 



Similarly, we can get the renormalization constants for 
the other octet baryons [21 
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I + 3/3eBi 
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C. Coefficients Crs, G,;b and Gk.b 

The coefficients G^s and GzB, where z is either e or 
K, were calculated in Ref. Q]. To express the results, it is 
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where e'^ — e^ = and e\ = 2, for Q^ = 0. The results 
for the coefficients are listed in Table IIIII (See Ref. [2| 
for details.) 



Octet baryon electromagnetic form factors witii 
pion cloud dressing 



D. 



From the current p9l) . using the expressions p9)) 
and ([20I) , and the coefficients in Table IIIH we can write 
down the Dirac form factors F\b for the octet baryons 
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^A |eoA + /3A(eos+ + eos" + cqe- )C'i 
+/5a(kos+ + Kos" + KoE- )-Di I , (38) 

^s{goE+ + (/?s+/3A)Si 

+ [^s(eos+ + eoE") + ^aSoa] C^i 

+ [^s(«os+ +«os") +^akoa]-Di|, (39) 

Zs jeos" 

+ [^s(eoE+ + eoE- ) + ^acoa] Ci 

+ [/3s(kos++koe-)+/5akoa]£'i}, (40) 

^s{SoE- -Ws+MBi 

+ [^s(eoE" + SoE- ) + ^acoa] Ci 

+ [/Ss(kos'' + KoE- ) + Pakoa] Di } , (41) 

Zs {SoHO + 2/3sBi 
+/3s{eos»+2eos-)Ci 

+/3h(koho+2koh-)^i}, (42) 

Zs {goH- - 2/?H^l 

+/3s{2ios» + eos-)Ci 
+/3h(2«oh"+«oh-)^i}. (43) 



F2p — Zn I Kop + 2-B2 

+ ieop + 2eon)C2 + (iiop + 2ko«)-D2 } , (44) 

F2n = Zn I KOn — 2^2 

+ (2gop + gon)C'2 + (2kop + iion)D2 } , (45) 



F2A = 



Za |koa + 13 A{eoT.+ + eos" + eos-)C'2 



(46) 



(47) 



-F2E0 



-P2E- 



+/5a('«oe+ + Kos" + Kos- )-D2 1 , 

i^2E+ = ^s {koe+ + (;9s + I3a)B2 

+ [/^E(eoE+ + eos") + /^acoa] C'2 
+ [/3e(kos+ + KoE") + /^akoa] -D2I 

■^E I '«0EO 

+ [/3E(eoE+ + eoE-) + /Jagoa] C2 

+ [/3e(koe++Soe-)+/3akoa]£'2}, (48) 

^e{^oe- -(/3e+/3a)S2 

+ [/3E(eoEO + eoE-) + /^agoa] C'2 

+ [/3e (koeo + KoE- ) + /3akoa] D2 I 

F2H0 = Zs{^oh«+2/3hS2 

+f3s{eoE»+2eoE-)C2 

+/3h(koho+2koh-)£'2}, 
F2H- = Ze{^oh--2/3hB2 

+f3s{2eoEO + eos-)C2 

+f3s{2kos« + iioE-)D2} ■ 



(49) 



(50) 



(51) 



The magnetic moment of the octet baryon member B 
is defined in terms of the magnetic form factor Gmb — 
Fib + F2B at Q"^ = 0, according to hb = Gmb{0)2M^- 
The results for the magnetic moments are usually ex- 
pressed in terms of the nuclear magneton fiN = „Jl , 

namely, ^s = G'Afs(0)^/iiv- 



Similarly, for the Pauli form factors F2B, we can write 



V. PARAMETRIZATIONS 

In the previous sections we have defined the general 
structure of the valence quark part and the functions for 
the pion cloud. We still need to specify the quark currents 
in terms of the functions fin{Q^) (i = 1,2, n = 0, ±), 



scalar wave functions iI>b{P, k), and the functions for the 
pion cloud effects, Bi, Ci, and Di {i = 1,2). 

A. Parametrization of the quark current 



To parametrize the quark current p4)) . we adopt the 
structure inspired by the VMD mechanism as in Refs. p], 
i, 



/lo = Aq + (1 - Ag) 2 I n2 + ^° 



MlQ^ 



/2± = K± i rf± 



.2 



t2 _i_n2 



(l-rf±) 



Ml 



Ml + Q2 [ ' 



^^"='^M'°^^^+('-'°)wTVr 



(52) 



where m„ , m^ and M^ are the masses respectively cor- 
responding to the light vector meson m„ ~ rup, the (j) 
meson (associated with an ss state), and an effective 
heavy meson with mass Mh = ^Mm to represent the 
short range phenomenology. For the isoscalar compo- 
nent it should be rriy = rrii^, but we neglect the small 
mass difference between the p and w mesons, and use nip. 
The coefficients Cq, c± and do, d± were determined in the 
previous studies for nucleon (model II) [I| and n~ Q- 
The values are, respectively, c+ = 4.160, c_ = 1.160, 
d+ = d- ^ -0.686, Co = 4.427 and do = -1.860 [|. The 
constant Xq = 1.22 represents the quark density number 
in deep inelastic scattering |l|. 

The isovector k+ = 2k„ — Hd and isoscalar k_ — 
-^{2ku + Kd) nucleon anomalous magnetic moments were 
adjusted differently in the past based on the other stud- 
ies: nucleon elastic form factors [1] or octet baryon mag- 
netic moments [2]. In Ref. [i| the u and d quark anoma- 
lous magnetic moments were fixed to describe the nu- 
cleon (proton and neutron) magnetic moments using the 
models without the contributions from the pion cloud at 
Q^ = 0. More recently, these parameters were updated 
in a model for the octet baryons with the pion cloud to 
reproduce the octet baryon magnetic moments [2| . That 
model was not constrained by the finite Q^ data, but 
only by the Q^ = data. In the present study because 
we want to describe also the finite Q^ data, we need to 
relax the conditions to fix k„ and Kd. The quark anoma- 
lous magnetic moments will then be constrained by the 
physical and lattice data. As for Kg (the strange quark 
anomalous magnetic moment) it was fixed as Kg = 1.462, 
to reproduce the f2~ magnetic moment [Z\. 

The quark form factors parametrized by a VMD mech- 
anism in Eq. ([5^ are particularly convenient to extend 
the model to the lattice QCD regime, because they are 
written in terms of the vector meson and nucleon masses. 
Then, the extension can be done replacing these masses 
by those of the lattice regime. 



B. Scalar wave functions 

Using the expressions P^ and (PU]) . we can determine 
all the electromagnetic form factors of the octet baryons. 
To do so, we need to specify the scalar wave functions "05 ■ 
Generalizing the form used previously for the nucleon [1| , 
we assume the scalar wave functions for the octet baryons 
in the form: 



MP,k) - 



Nn 



moifii +Xjv)(/32 +Xjv)' 

Nk 

moilii +Xa)(/33 +Xa)' 



(53) 
(54) 
(55) 



^H(P,fc) = ,. ^ ~. ^ . , (56) 



where Nb {B = N, A, S, .=,) are the normalization con- 
stants, and 



(Mb - mpf - (P - kf 
Mbtud 



(57) 



Note that, except for the masses, the A and S scalar wave 
functions are the same. The normalization constants Nb 
are determined by 



|^s(P,A:)|2 = l, 



(58) 



where P — [Mb, 0, 0, 0) is the baryon four-momentum at 
its rest frame. 

In Eqs. ([55|) - ([55)) the parameters (ii (i = 1, .., 4) define 
the momentum range in units of mu.^ In the scalar wave 
functions ([551) - ([5S)) with the assumption /32,/33,/34 > /3i, 
we associate /3i with the long range scale (low momentum 
range) that is common to all the octet baryon members, 
and the remaining Pi with the shorter spatial range scale. 
The parameters /3i and /?2 can be determined by the nu- 
cleon data only (no strange quarks) in a model without 
the pion cloud. (See e.g., Ref. fl|.) The parameters P^ 
and Pi are associated with the strange quark. While P^ 
is related to the system with one strange quark and P^ 
is related to that of the two strange quarks. Since the 
strange quark is heavier than the u and d quarks, we can 
expect p2 > Ps > P4- The parameters pi {i = 1, .., 4) will 
be determined later. 



^ In the baryon B rest frame, it reduces to 

ft + Xs =(ft-2) + 2-^. 



C. Extension of the model for the lattice regime 



parametrization: 



We now discuss the extension of the model for the lat- 
tice regime. In this regime the mass of the octet baryon 
B is characterized by the pion mass in lattice mj^''** and 
denoted by M^**. The current j'^ [Eq. P^ ] is also char- 
acterized by the corresponding pion mass in terms of the 
two components ji and J2, which are represented based 
on the VMD parametrization in Eq. ([5^ . 

In the quark current (J14l) we replace the coefficient 
of the Pauli form factor 1/{2Mn) by 1/(2M^") in the 
lattice regime. As for the quark form factors, we use 
Eq. ([5^ with the meson masses replaced by the respec- 
tive lattice masses as in the previous studies [3, |20, IsJ ■ 
Namely, we replace my by the lattice p mass m'"*', and 
the effective heavy meson mass of Mh = 2Mm by 2M^**. 
On the other hand, the physical mass is used for tti^, 
since presently lattice simulations are performed using 
the physical strange quark mass. To represent the p me- 
son mass in the lattice regime, we use the following ex- 
pression based on the lattice studies made in Ref . l49| , 



Si = Bi 



^latt _ „ I „ (flatty 



(59) 



where ao = 0.766 GeV and aa = 0.427 GeV-^. With 
this procedure we can define unambiguously the quark 
current. 

As for the wave function ^s, AIb is replaced by M^** 
in the lattice regime. This applies for the scalar wave 
functions ([551) - ([5S)) . There is no need of modifying the 
diquark mass in the lattice regime, since the electromag- 
netic form factors are independent of it. We assume that 
the range parameters Pi {i = 1,..,4) are independent of 
the baryon masses and therefore independent of the lat- 
tice pion mass mj^**. We can expect this approximation 
to work for a certain range of the mj^"** values, but it 
breaks down for the larger values of mj^"'* [S^ ■ 

Using the model extended to the lattice regime, 
namely, using the quark currents and baryon wave func- 
tions in the lattice regime, we can calculate the form 
factors Fib and F2B in the lattice regime via Eqs. P^ 
and (|20p . using the lattice regime masses corresponding 
to mJr"**. However, note that the results include only 
the valence quark contributions, but not the pion cloud 
contributions. 



D. Pion cloud factors 



To describe the pion cloud contributions for the octet 
baryon electromagnetic form factors, we use the following 






Ci=Bi 



A?2 '' 



A?2 + Q^ 



^1 = ^17X2 
B2 = B2 



, Q^K\, 



(A?3 + Q')^ 

Ail 



C'2 — Co 



Aii + g^ 

Aja 
Ai2 + Q^ 



A2 



^^-^^',Al3 + Q^ 



(60) 
(61) 
(62) 

(63) 
(64) 
(65) 



The parametrization is phenomenological and motivated 
by the expected falloff of the quark-antiquark contri- 
butions at very large Q^ [50], as well as the magni- 
tude of the pion cloud contributions estimated for the 
77V — >• A reaction |18l - [2(j | . The pion cloud contributions 
should falloff by a factor l/Q^ faster than the falloff of 
the valence quark contributions. In principle, Bi and 
B2 should contain the pion electromagnetic form factor 

FjriQ^) ~ ( 1 + 05 ) ' ^^* ^^ adopt simplified functions 
using just one cutoff parameter with the powers 4 and 5, 
respectively. 

In the above Bi (= Ci), B2, C2 and D2 represent the 
values of the respective functions at Q^ =0. The values 
were determined previously by the octet baryon magnetic 
moments in Ref. ^. D[ is a new constant defined as 
D[ = -^^k{0). However, we do not use the values 

determined in Ref. [3] in this study. Instead, we will 
determine them directly by the nucleon form factor data 
at finite Q^, and by the A, S and S physical magnetic 
moment data, after determining the contributions from 
the valence quarks. The adjustable parameters of our 
pion cloud parametrization are then, the coefficients for 
Q^ — 0, and the cutoffs Aii and Aai (i — 1, 2, 3). But for 
simplicity, we use and vary only two independent cutoff 
values for the Dirac and Pauli form factors, respectively 
Ai and A2 in this study. 



E. Separating the pion cloud contributions 

Based on the discussions in the previous sections, the 
baryon form factors Fis {i = 1,2) may be decomposed 
into 

Fb{Q^) = Zb [Fob{Q^) + SFwiQ^)] , (66) 

where Fiqb = eos, F20B = kob and corresponding 
SFi^2B the pion cloud contributions, which are given in 
Eqs. (p6l) -(|5T |) . We may regard that ZbFmb represent 
the effects of the valence quarks, and ZbSFib those of 
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the pion cloud. The same decomposition can be apphed 
for the electric and magnetic form factors: 



Geb[Q'')^Zb[Geob[Q'')^ 
Gmb{Q ) = Zb [Ga/osIQ ) 



5Geb{Q'')\ , 
h<5GMB(Q2)],(67) 



where Geob ^ eos-rKos, Gmob = eos + Kos, 5Geb = 
SFiB — t5F2b, and SGmb = SFib + SF2b- In this 
case ZbSGeb, and Zb&Gmb reflect the dressing of the 
pion cloud. To estimate the pion cloud contributions, 
we will compare the full result, Geb or Gmb, with the 
total contributions of the valence quark core, ZbGeob 
or ZbGmob- The difference is the pion cloud contribu- 
tions, ZbSGeb or ZbSGmb- Hereafter, we will use Gxb 
to express the electric form factor Geb {X ~ E) or the 
magnetic form factor Gmb {X = M). 

Note that we can alternatively define the contributions 
from the pion cloud as the difference between the bare 
form factor, Geob or Gmqbi and the full form factor, 
Geb or Gmb, instead oi ZbGeob or ZbGmob described 
above. We will refer to these terms as the effective pion 
cloud contributions. By this definition we may have small 
effective pion cloud contributions in the cases, e.g., where 
ZbSGeb or ZbSGmb are significant^. Such an example 
is the model in Ref. [3|. There the bare contributions 
for the nucleon magnetic form factor at Q^ = 0, given 
by ZnGmon, was about 60% (namely, the pion cloud 
contributions were ~ 40%), although the effective pion 
cloud contributions were only 5%. 



VI. RESULTS 

In this section we determine the parameters of the 
model, and then present results for the valence quark 
contributions for the octet baryon electromagnetic form 
factors in the lattice regime, as well as those for the phys- 
ical regime which includes the pion cloud effects. While 
the expressions related to the valence quark contribu- 
tions have been presented in Sec. IIII Dl the formalism 
related to the pion cloud dressing has been presented in 
Sec. UVDl 

The parameters associated with the valence quark de- 
grees of freedom and the ones associated with the pion 
cloud dressing, are determined by a global fit to the octet 
baryon electromagnetic form factors, namely, the lattice 
data, nucleon physical data, and the octet baryon physi- 
cal magnetic moments. Details are described next. Once 
we have fixed the relevant parameters, we will discuss the 
results for the valence quark contributions and the effects 
of the pion cloud dressing. 



^ The effective pion cloud contributions for the form factor Gxb 
is ZgSGxB ~ (1 ~ Zb)Gxob- If ^B differs substantially from 
1, the pion cloud contributions ZbSGxb will be modified by the 
term —(1 — Zg)GxoB, and large cancellation can happen. 



A. Global fit 

In a regime where the pion cloud effects are small, such 
as the lattice QCD simulations with heavy pions, or the 
physical regime at high Q^ , the electromagnetic form fac- 
tors of the octet baryons should be well described only by 
the valence quark degrees of freedom. Then, we can use 
the lattice QCD data with large pion masses to calibrate 
our model extended to the lattice regime based on the 
formalism described in Sec. IV CI without the pion cloud 
effects. 

To describe the physical octet baryon systems, only 
the valence quark degrees of freedom are usually insuffi- 
cient, and explicit pion cloud effects are necessary. Ex- 
cept for the magnetic moments [4] , there are no physical 
data available for the octet baryon form factors besides 
the nucleon system. Therefore, the nucleon system in 
the physical regime, is ideal to study the pion cloud ef- 
fects on the form factors at finite Q^. In this case the 
nucleon form factors can be described by the mixture of 
the valence quark and pion cloud contributions given by 
Eqs. dSH), (EZl), dH and gSl). Also the A, I]± and S^'" 
magnetic moments, defined by Eqs. ([38l) - (|43)) and ([46]) - 
([STj) at Q^ = 0, can be used to constrain the effects of 
the pion cloud at Q^ = 0. 

Summarizing, to adjust the parameters of our model, 
we perform a fit to the lattice data for the valence quark 
part by extending the model to the lattice regime, and 
for the pion cloud contributions we perform a fit to the 
physical data of nucleon form factors and octet baryon 
magnetic moments. In the latter, the physical regime, the 
expressions for the form factors are given by extrapolat- 
ing the bare form factors to the physical case (m^ = 138 
MeV) , and the pion cloud contributions are parametrized 
according to Eqs. (f60 |) -([65 l) . Next, we describe specifi- 
cally the lattice QCD data and the nucleon physical data. 
As for the magnetic moments, we use the experimental 
data for the A, S^, and S°^~ [5l[. We do not include the 
neutron and proton magnetic moments in the fit, since 
the nucleon magnetic form factor data are also included 
in our analysis. 



1. Lattice data 

The lattice data adopted in this work, are the octet 
baryon electromagnetic form factor data from Ref. [52l |. 
where n, p, S*, 'EP'~ form factors were calculated sys- 
tematically at finite Q^ for the first time. The data from 
Ref. [52] are composed of four sets of unquenched simu- 
lations associated with the pion masses 354, 495, 591 and 
680 MeV, but restricted to the region Q^ < 1.5 GeV^. As 
a total we have 136 data points for both Geb and Gmb- 
To extend our model (valence quark contributions) to 
the lattice regime, we follow the procedure described in 
Sec. lVCl The relevant variables necessary are the masses 
associated with the lattice QCD simulations. The corre- 
sponding values are presented in Table IIVI according to 
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FIG. 3: Bare electromagnetic form factors for the proton determined by the global fit, compared with the lattice data from Ref. |52|| . 
The lines are the lattice regime (solid line) and the physical regime (dashed line). For m,r = 138 MeV, we include the physical data. The 
thin solid lines in the bottom panels do not include the pion cloud contributions. 



the simulations of Ref. [52 
mined through Eq. ([55)). 



and the 



values deter- 



2. Nucleon physical data 

For the nucleon we have included the proton electric 
and magnetic form factors {Gep and Gmp)-, and neu- 
tron electric and magnetic form factors {Geu and Gmti)- 
Since one of our goals in this work is to describe the octet 
baryon lattice data for Geb and Gmb, we do not adopt 
the ratio Geb/Gmb, although it is considered in many 
studies of the proton electromagnetic form factors. 



The proton data can be extracted using the classical 
Rosenbluth separation technique from the cross section 
data, or the polarization transfer method developed at 
Jefferson Lab to measure the ratio Gep/Gmp [53, [5J]. 
The results of the two methods show discrepancies that 
can possibly be explained by including the two-photon 
exchange corrections in the results of the Rosenbluth 
separation method [iTJ, l55| . An analysis that takes into 
account the two-photon exchange corrections and uses 
the cross section information to determine the values of 
Gep and Gmt, sep arately (not just the ratio), was pre- 
sented in Ref. [55|. Since our calculations are performed 
in the impulse approximation, we compare our results 
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FIG. 4: Bare electromagnetic form factors for the neutron determined by tlie global fit, compared with the lattice data from Ref. [52l |. 
The lines are the lattice regime (solid line) and the physical regime (dashed line). For m^ = 138 MeV, we include the physical data. The 
thin solid lines in the bottom panels do not include the pion cloud contributions. 



with the analysis of Arrington et al. [55], where the two- 
photon exchange contributions were subtracted. To in- 
clude the recent, high Q^ results for Gep/Gmp from Jef- 
ferson Lab [5J], we convert the ratio Gep/Gmp into Gep 
using the fit to the Gj\/p presented in Ref. [S^ • The fit is 
accurate for a large Q^-range {Q^ =0—10 GeV^). Over- 
all, we have 50 data points for Gep and 56 data points 
for Gmp- 

As for the neutron form factors, we collect the data 
from different groups. We prefer to use the data ex- 
tracted from a deuterium target rather than those ex- 
tracted from a "^He target, since the former is expected 



to have fewer nuclear corrections. For Geu we use the 
data from Mainz ^, NIKHEF ^, MIT-Bates [58^ and 
Jefferson Lab _59, 60]. The results from Ref . i60] obtained 
using the '^He target, corresponding to the highest Q^ re- 
sult for Geti (Q^ = 3.4 GeV^), are also included. Also the 
results obtained from the deuteron electric quadrupole 
moment are included 1611 . For Gmu we adopt the data 
used by Bosted et al. ^ in the global fit of the nucleon 
data, as well as more recent data from Mainz [63| and 
Jefferson Lab ]64]. Totally, we have 29 data points for 
Geu and 67 data points for Gmu- 

The inclusion of the high Q^ region nucleon data. 
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FIG. 5: A bare electromagnetic form factors determined by the global fit. The lines are the lattice regime (solid line) and the physical 
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m^(GeV) 


mp(GeV) Miv(GeV) ME(GeV) 


MH(GeV) 


0.138 


0.779 


0.939 


1.192 


1.318 


0.351 


0.820 


1.150 


1.349 


1.438 


0.495 


0.871 


1.290 


1.410 


1.475 


0.591 


0.915 


1.366 


1.448 


1.491 


0.690 


0.964 


1.490 


1.524 


1.546 



TABLE IV: Masses of the octet baryons {M^ , My. and Ms) 
obtained in lattice QCD in Ref. [52 ]. The p meson mass is 
obtained using the parametrization (I59p . from Ref. [43|. The 
first row corresponds to the physical values. In the physical 
case it is also A/a = 1.116 GeV. 



where pion cloud effects are expected to be small, is 
important for the cahbration of our model in the va- 
lence quark sector, since the information of lattice sim- 
ulations is restricted to the low Q^ region. We can test 
the parametrization directly by the nucleon elastic form 
factor data that are extended for the proton up to 9 and 
31 GeV^ respectively for Gep and GMp^ and for the neu- 
tron up to 3.4 and 10 GcV^, respectively, for Geti and 
Gj\/„. 



3. Details of the fit 

The parameters associated with the valence quark con- 
tributions are adjusted to the lattice data (bare form 
factors), the bare part of the nucleon form factors at 
the physical point, and octet baryon magnetic moments. 
They are the momentum range parameters of the scalar 
wave functions, /?i {i = 1, ..,4), and the u and d quark 
anomalous magnetic moments k„ and Hd, respectively. 
The strange quark anomalous magnetic moment Ks is 
kept unchanged, since it was fixed by the 57^ magnetic 
moment, which is not affected by the pion cloud (3|. 

The pion cloud contributions are adjusted by the nu- 
cleon physical data and the octet baryon magnetic mo- 
ment data. We start by calculating the bare form factors 
cqb and kqb at the physical point (rriTr = 138 MeV), us- 
ing only the parametrization for the bare form factors. 
Next, we add for each form factor, Fi and F2, or alterna- 
tively Ge and Gm , the pion cloud contributions given by 
Eqs. ([Sn)) - (|55|) . For this, we need to fix the parameters 
Bi,B2,Ci,D[,D2 and the cutoffs Ai and A2. These pa- 
rameters are adjusted by the nucleon form factor data for 
Ge and Gm, and octet baryon magnetic moments, since 
the magnetic moments are proportional to Gmb{0)- In 
the latter we have additional constraints for the coeffi- 
cients Bi,Ci, C2 and D2 [2| • (The calibration of D[ , 
Ai and A2 can be done only with Q^ > data). Once 
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FIG. 6: S+ bare electromagnetic form factors determined by the global fit. The lines are the lattice regime (solid line) and the physical 
regime (dashed line). The lattice data are from Ref. |52ll . 



we have fixed the parameters associated with the pion 
cloud, we can obtain explicit expressions for the pion 
cloud contributions for the octet baryon electromagnetic 
form factors. 



For very large pion masses, it is not reasonable to as- 
sume the baryon wave functions can be described by the 
same parameters, but some dependence on the baryon 
masses (indirect dependence on the quark mass) may en- 
ter. Also, for small pion masses (ttItt < 400 MeV) one 
can expect that lattice simulations may be affected by 
the pion cloud effects, and as a consequence form fac- 
tor data should differ from the calculation based solely 
on the valence quark degrees of freedom. The ideal sit- 
uation would be to use several sets of lattice data in the 
region 400 MeV < rn^ < 600 MeV. In this way we expect 
to avoid both the pion cloud contamination, and the very 
large rriTr region where the model can fail. In the previous 
studies, this way of extension to the lattice regime was 
very successful [20, [sj . 



A preliminary analysis of the lattice data has revealed 
that the set tti^ = 591 MeV is particularly difficult to de- 
scribe in our model (large x^ per data point as presented 
in Table |V| , although the quality of the fit is improved 



for the next set"^ {m.,^ = 691 MeV). Using efficiently the 
available lattice QCD data for the octet baryon electro- 
magnetic form factors, and at the same time to keep a 
reasonable description of the present model, we perform 
a fit to the first two sets of lattice data, m^r =351 and 
495 MeV. (The risk of the pion cloud contamination in 
the m^ = 351 MeV set is compensated by the increase of 
accuracy in the set.) 

Another point to be noted in our analysis is the lat- 
tice data for the neutral particles, n and "EP . The results 
for the n and S" electric charge form factors at Q^ = 
(their charges) differ from zero, and this fact suggests 
that the lattice results have some systematic errors. This 
may be a consequence of incomplete cancellations among 
the contributions from the different quark flavors for the 
form factors [52l |. and we have to take this into account. 
In order to use the n and 'EP data but achieve a reason- 
able accuracy, we reduce the respective impact on x^ by 



^ Although the quality of a global fit is measured by the total 
X^ divided by, [(the number of data points) - (the number of 
parameters)] considered in the fit, to estimate the quality of the 
subsets of data, it is simpler to ignore the subtraction of the 
number of parameters, since the parameters are not fixed by one 
subset of the data. 



15 




1 

0.8 

"^ 0.6 

a 

O 0.4 

0.2 





s\ m =495MeV 



2-1 • Lattice: m = 464 MeV 

J , , , , I , , , , \ , , , , \ ^ 

0.5 1 1.5 



- 






m =354 MeV 

7t 



-0.04 



• Lattice: m = 372 MeV 



0.5 1 1.5 

Q' (GeV') 



1 

0.8 

r^^ 0.6 

O 

O 0.4 

0.2 





1 1 1 1 1 

V 

s 

: e" 

1 1 1 1 1 


m =354 MeV 

71 


• Lattice: m^= 372 MeV 


" 



0.5 1 1.5 



FIG. 7: E" bare electromagnetic form factors determined by the global fit. The lines arc the lattice regime (solid line) and the physical 
regime (dashed line). The lattice data are from Ref. [65| for m^ = 372 and 464 MeV. 



doubling the respective statistical errors. 

As already mentioned, the pion cloud parametrization 
is calibrated using the A, E='=, S"'" physical magnetic 
moment data and the nucleon physical data. Some of 
the octet magnetic moment data are extremely accurate 
with error bars of less than 1% (A and E~). A mini- 
mization of x^ with such small error bars of Gmb (0) will 
impose strong contributions on x^ in the Q^ = region, 
and reduce the relative impact of the region Q^ > 0, 
represented by the nucleon form factor data. To achieve 
a good description (low x^/^p; where Up is the number 
of the data points) for the nucleon data, we reduce the 
weight of the octet baryon magnetic moment data, by 
doubling the experimental errors. 

We define the best model as the model that minimizes 
the total x^ associated with the lattice data, the octet 
magnetic moment data and nucleon physical data as de- 
scribed before. Totally, we have 272 lattice data points, 
5 experimental magnetic moment data points, and 202 
nucleon physical data points. 

A fit with no constraints leads to a good description 
of the lattice data (small y^ per data point) , but a poor 
description of the nucleon physical data with x^ P^r data 
point > 2. That fit generates also a large contributions 
from the pion cloud effects, and has a large extension in 
Q^ from the pion cloud effects, in particular for Gm due 



to a large cutoff A2 compared to 1 GeV. We interpret 
this result as a consequence of the dimension of the lat- 
tice database (272 data points) for 6 baryons including 
the nucleon, to be compared with the nucleon physical 
data (202 data points). This procedure reduces the rel- 
ative importance of the nucleon physical data. Since we 
want to describe well the nucleon system but simultane- 
ously to have a reasonable description of the lattice data, 
we reinforce the impact of the nucleon physical data by 
doubling the contributions of their x^ in the global eval- 
uation of x^ . The fit with this constraint leads a good 
qualitative description of the nucleon physical data (x^ 
per data point ~ 1-9), and also reasonable values for the 
cutoffs Ai and A2, where the values are smaller than 1 
GeV, or closer to 1 GeV, consistent with the pion cloud 
effects restricted to low Q^ region. The values of the 
cutoffs will be discussed in more detail later. As we will 
see, the final fit is consistent with the small pion cloud 
contributions for the octet baryon form factors. 

The quality of the fit, measured in terms of x^ per 
data point, can be understood from Table |Vl In the 
table we represent for each set, x^{Geb) in the first row, 
X^{Gmb) in the second row, and the combined result in 
the last row. With bold face ==l In boldface In boldface 
we represent the x^ for the sets considered in the fit. 
The rows associated with the physical regime (tott = 138 
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FIG. 8: S bare electromagnetic form factors determined by the global fit. The lines are the lattice regime (solid line) and the physical 
regime (dashed line). The lattice data are from Ref. [52l| . 



MeV) are the only ones that reflect the effects of the pion 
cloud. Note that the sets m^ = 591 and 680 MeV, are not 
included in our fit. In the last column we represent the 
partial ^{Ge) and 'yd^Gu)-, and the total x^j associated 
with respective sets. In this case the contributions from 
n and "E? are not included. 



B. Bare octet form factors 

The results of the fit for the octet baryon electromag- 
netic form factors are presented in Figs.l3llT0l The values 
of the pion mass in the fit are mj; — 354 and 495 MeV. 
For the nucleon [Figs. [3] and |3], we present also the phys- 
ical case (to^ = 138 MeV). The parameters associated 
with the results are 



/3i = 0.0440, /32 = 0.9077, 
iiz = 0.7634, /34 = 0.4993, 
K„ = 1.6690, Kd = 1.9287. 



(68) 



We note that the (bare) quark anomalous magnetic mo- 
ments, Ku and Kd, are similar to the previous results 
obtained in Refs. [l|, Q], within a 16% variation. (In 
Refs. [l|, 0] the values obtained are k„ = 1.778 and 
Kd = 1.915, while in Ref. Q they are Ku — 1.929 and 
Kd = 1.919.) 



Furthermore, the values for /3i and /32 are also sim- 
ilar to those of the model in Ref. |l| (/3i — 0.049 and 
/?2 — 0.717). Interpreting /3i as the spatial long range pa- 
rameter common to all the systems associated with one 
light quark (TV, S, A and S), P2, P3 and /34 are interpreted 
as spatial short range parameters associated with spatial 
extension of the qiqi, qis and ss quark pairs respectively, 
where qi stands for a light quark. In this respect we 
should expect /32 > /^a > /?4, consistently with Eq. (68), 
a decreasing of the spatial extension of the systems, grad- 
ually from the nucleon followed by the A, T, and finally 
by the S system. The systems with one strange quark 
are more compact than the ones with no strange quarks, 
and the systems with two strange quarks are even more 
compact than the ones with only one strange quark. We 
will return to this issue later when we discuss the charge 
and magnetic squared radii. 

In Figs.[51fTUlwe compare the results of our fit for the n, 
p, A, S°'='=, and S°'~ with the lattice data for m^ = 354 
and 495 MeV. The lattice data [53] used in the fit are 
shown in Figs. [H IH [6l [U [H andfTOl respectively, for the 
p, n, S"*", E~, 'EP, and S~. For the proton and neutron 
we include also the physical case {itIt^ = 138 MeV) and 
compare the results with the experimental data (details 
will be discussed later). We have no data for the A and 
E° from Ref. [53], respectively, shown in Figs. and [3 
but we compare our results with the quenched lattice 
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FIG. 9: S bare electromagnetic form factors determined by the global fit. The lines are for the lattice regime (solid line), and the 
physical regime (dashed line). The lattice data are from Ref. |52|| . 



QCD simulation results from Ref. [65J with the closest 
pion mass values, respectively m^r = 372 and 464 MeV, 
for the single point at Q^ = 0.23 GeV^ calculated in that 
work. 

Overall, we have a very good description of the data 
for p, S and even S~. As mentioned already, the re- 
sults for n and "EP should be taken with caution, since 
the lattice simulations have systematic deviations from 
the expected results, particularly for the electric charge 
form factors. Nevertheless, we have a good global de- 
scription of the nucleon and E systems. As for the S, 
the lattice results are closer to the extrapolation to the 
physical limit (dashed line) than the calculation in the 
lattice regime (solid line). This can be caused by the 
poor quality of the Ge data for S° as discussed before, 
or a limitation of our simplified approach. Further lat- 
tice QCD simulation data, consistent with the n and 'EP 
charges, are necessary to clarify this point. As for the 
other neutral particles, A and E", our predictions are 
consistent with the simulations of Boinepalli et al. [6a|, 
within about 2 standard deviations, for the closest pion 
mass values used. A final remark is on the difficulty of 
the present approach in describing the n and S° data. 
Since we try to describe well Fi and F2 (or Ge and Gm) 
simultaneously, inaccuracy in one set (say Ge) will affect 
the description of the other set (say Gm)- 

For the nucleon system we compare also the results 



extrapolated to the physical limit of the bare form fac- 
tors with the physical data (see the bottom panels with 
TOtt = 138 MeV in Figs. E] and H]). Note that these results 
should not be compared directly, since we have not yet 
included the pion cloud effects in the calculation. From 
the figures, our aim is to see whether or not the pion 
cloud effects are indeed important, and if reasonable de- 
scription of the data can be achieved without the pion 
cloud effects. We plot then the results of two different 
calculations. The first one is the extrapolation to the 
physical limit of the model, represented by the dashed 
line (the same as in the upper panels). The result is ob- 
tained setting Zjsi = 1, in the expression for the nucleon 
form factors [Eqs. dMl, dSZl), &, and (gS])], and remov- 
ing the pion cloud contributions. Plotted in the same 
figures (thin solid line), is the same calculation but for 
the case Z^ =0.885, given by the fit, the contributions 
exclusively from the valence quark core. From the results 
shown in the figures, we conclude that the data are well 
described by a model with no pion cloud effects, although 
in the region of high Q^ (say Q^ > 1 GeV^), the model 
with Zn = 0.885 given by the thin solid lines (solely from 
the valence quark contributions), are closer to the data 
than those of Zm = 1 • 
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FIG. 10: H bare electromagnetic form factors determined by the global fit. The lines are the lattice regime (solid line), and the physical 
regime (dashed line). The lattice data are from Ref. [5^ . 



C. Pion cloud contributions 

We discuss now the calibration of the pion cloud ef- 
fects. The values in the paranietrization associated with 
the pion cloud effects are 

Bi = 0.04343, B2 = 0.21477, C2 = 0.02266, 
D[ = -0.17637, L»2 = 0.08551, (69) 

with the cutoff values, 



is the neutron electric form factor. In the figure, the 
contributions of the pion cloud are represented by the 
bands that fill the difference between the valence quark 
contributions {ZbGxqb) and the full result (Gxbi solid 
line) . 

Observing the pion cloud contributions for the nucleon 
electromagnetic form factors in Fig. [Til '^e conclude that 
the contributions are similar for both the proton and neu- 
tron magnetic form factors. In both cases contributions 
amount to 10-14% in the region of Q^ = - 0.5 GeV^, 
and fall to less than 5% around Q^ = 2 GeV^, and even 



Ai = 0.7732 GeV, A2 = 1.2455 GeV. (70) become less than 1% for Q^ > 5 GeV 



The quality of the fit associated with the pion cloud ef- 
fects is measured by the partial x^ values for the nucleon 
system at the physical point, given by the last column 
in Table FVl (x^ per data point — 1.93). Then, we can 
conclude that the nucleon data are described better than 
the lattice data (x^ per data point of 2.9 and 5.2 for the 
sets TO^ = 354 and 495 MeV, respectively). 

Since we cannot isolate the pion cloud contributions 
from the valence quark contributions in the experimental 
data, we analyze the pion cloud effects by comparing the 
individual components of the nucleon form factors with 
the full result. The results for the nucleon are presented 
in Fig. llli where the form factors are renormalized by the 

dipole form factor Go = ( 1 + ^fj 1 . The exception 



The analysis for the electric form factors is more deli- 
cate. For the proton there are w 12% contributions from 
the pion cloud near Q^ = 0, and they fall to 1% near 
Q^ = 1 GeV^, and stabilize to 5% negative contributions 
for Q^ « 5 GeV^. In the larger Q^ region one must be 
careful, since Ge approaches zero and the ratio is not 
meaningful. For Q^ = 10 GeV^ the valence quark con- 
tributions are larger than 90%. As for the neutron near 
Q^ = 0, where GEn{0) = 0, the pion cloud contributions 
dominate. Near Q^ — 1 GeV'^, the pion cloud effects are 
about 10% and drop to less than 4% for Q^ = 4 GeV^, 
and even smaller for larger Q^. For Q^ = 10 GeV^ the 
valence quark contributions dominate to give more than 
98%. 

The slow falloff of the pion cloud contributions for the 
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1.600 
1.857 
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2.273 
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1.700 
2.083 
1.933 



TABLE V: x^ decomposition of the best fit to the lattice data and nucleon physical data. The numbers in the rows corresponding 
to m-n — 138 MeV (nucleon data) include the effects of the pion cloud. See also the discussion in the text. For the lattice data 
X^ per data point is 5.00 (2.93 excluding the n and EP data). 



electric form factors compared with those for the mag- 
netic ones, is due to the enhancement of the F2 contribu- 
tions for Ge by the pre- factor Q^, and the function form 
for the pion cloud contributions. Since the pion cloud 
contributions are regulated by the cutoff A2 — 1.24 GeV 
(Al = 1.55 GeV2) which is larger than Ai ~ 0.77 GeV 
(Af = 0.59 GeV^), the electric form factor is extended to 
the higher Q^ region due to the range of the pion cloud 
contributions in F2. 

We point out that the qualitative differences of the 
pion cloud contributions between the present results and 
those of the previous study in the octet baryon magnetic 
moments 2]. In the previous study the Q^ dependence 
was not taken into account (no use of lattice data, nei- 
ther the nucleon physical form factor data), and the pion 
cloud contributions were significantly larger. Thus, we 
can conclude that the Q^ dependence of the pion cloud 
contributions is very important to constrain the pion 
cloud contributions, in particular for the nucleon system. 



namical behavior (Q^ dependence) of the form factors, we 
include also the lattice QCD simulation data correspond- 
ing to the lowest pion mass from Ref. [6a| (TOtt = 306 
MeV). In principle, the lattice data should be compared 
with the bare form factors (dashed line), unless strong 
quenched effects are expected. 



An important conclusion from the figures is that the 
magnitude of the pion cloud contributions is small. Based 
on the magnetic moments, we have more significant con- 
tributions from the pion cloud for E~, E"*", and A, respec- 
tively about 21%, 14% and 14%, and less than 10% for all 
other cases. With the exception of the neutral particles 
A, T,^, and EP, where pion cloud effects can be dominant, 
the pion cloud contributions for the electric form factors 
are smaller than those for the magnetic form factors. 



D. Dressed form factors 



Taking into account the results for the bare form fac- 
tors extracted from the lattice data, and the pion cloud 
parametrization from the previous section, we use the 
expressions in Sec. lIVDJ to predict the dressed, physical 
form factors for the A, E, and S systems in the physical 
regime (m^ = 138 MeV). The predicted results are pre- 
sented in Figs. [13 [T2] and [131 In the figures we show the 
full results (solid line) and bare results (dashed line) . In- 
cluded in the figures are also the magnetic moments when 
known (A, E*, and S°'^). To have an idea for the dy- 



Overall, our predictions for the octet baryon electro- 
magnetic form factors, as functions of Q^, are consistent 
with the results for the magnetic moments (at Q^ = 0), 
and close to the lattice QCD simulations. The major ex- 
ception is the results for S°, where we observe the clear 
deviation from the experimental result for Gm{0), and 
also from the lattice data. We recal again that this can 
be a consequence of the difficulty in describing the S lat- 
tice data. [See Figs. [51 and [TUl] 
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FIG. 11: Nucleoli electromagnetic form factors including the sum of the bare and the pion cloud (bands). Data are from Refs. [53h64| | . 
See text for details. 
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FIG. 12: A electromagnetic form factors for the total result (solid line) and the bare result (dashed line). The data point for Q^ = is 
the result of the magnetic moment J5l| . The lattice data point (filled triangle) for Q^ = 0.23 GeV^ is from Ref. [65l | (m-^ = 306 MeV). 



E. Electric charge and magnetic dipole radii 

The electric charge squared radius for a charged parti- 
cle is usually defined as** 



(-1) = - 



6 



dG 



EB 



Geb{0) dQ^ 



(71) 



Q2=0 



For a neutral particle the same expression can be used 
but setting Geb{0) -> 1- The definition ([7T|) has advan- 
tages for comparing the radii of particles with different 
charges such as p and E~, and one can relate the corre- 
sponding baryon electric charge radii. As for the mag- 



** Some authors 16511 exclude the factor Geb{0) from the {r^} 
definition. 
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FIG. 13: S electromagnetic form factors for the total (solid line) and the bare result (dashed line). The data point for Q^ = is the 
result of the magnetic moment [5J |. The lattice data point (filled triangle) for Q^ = 0.23 GeV^ is from Ref. [65l | {m-^ = 306 MeV). 



netic dipole squared radius, the most common definition^ 

is 



(rh) = 



6 



dG 



MB 



Gmb{0) dQ^ 



Q'^=0 



(72) 



We assume in this case that Gmb{0) is not zero, neither 
very small. The results for the electric charge squared 
radii and the magnetic squared radii are respectively 



^ Some authors |65t define {r'^j} without the factor G]v/b(0), but 



3(0) 



to compare the values of different baryons. 



presented in Tables IVII and IVIII (see columns {r"^) and 
(r^)). Experimental values |5ll. [GGhTII I are also included 
in Table IVTl for {r%), and in the caption of Table fVIII for 

Since in the present approach we can identify the va- 
lence quark (bare) contributions and the pion cloud con- 
tributions in the form factor Gxb {^ = E, Af ), we follow 
Eq. (|67|) and decompose Gxb into 



Gxb{Q^) 



G'xBiQ') + GlBiQ' 



(73) 



where G^xeiQ') - ZbGxob{Q') and G^s(Q') = 
ZbSGxbIQ^), are respectively the bare and pion cloud 
contributions. Based on the decomposition (j73l) and the 
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FIG. 14: S electromagnetic form factors for the total result (solid line) and the bare result (dashed line). The data point for Q = is 
the result of the magnetic moment [Sll . The lattice data point (filled triangle) for Q'^ = 0.23 GeV'^ is from Ref. |65|| (r7i,r = 306 MeV). 



definitions of radii (|7T|) and ((72)) , we can write 



where 



{rl\ 



-Zf 



6 



dG 



BOB 



Hi), = -Zi 



Gbb(O) dQ2 
6 cIGmob 



Q2=o 



Gmb{Q) dQ^ 



(74) 
(75) 

(76) 
(77) 



and (f|;) and (''^^ ) can be defined in a similar manner, 
but using 5Geb{Q^) and 5Gmb{Q^), or from Eqs. (|74)) 
and (|75p. by subtracting the bare components from the 
total. 

Since our form factors are determined numerically, we 
calculate the octet baryon electric charge squared radii 
and magnetic dipole squared radii, as well as the respec- 
tive components, using numerical derivatives. The con- 
tribution from each component is also presented in Ta- 
bles Ell and |VlTl 

From the global results for (t*!;) and (t\i\, we can con- 
clude that the nucleon system has larger spatial charge 
and magnetization distributions than that of the E sys- 
tem, and that the S system has the larger spatial charge 
and magnetization distributions than those of the S sys- 
tem. We will leave the neutral particles (n, A, YP and S°) 



out of the discussion for {r%). Our results suggest that 
the electric charge squared radii of 0.7-0.8 fm^ for pro- 
ton (p), 0.6-0.7 fm2 for S±, and 0.4 fm^ for S^. These 
results may support the general idea that the systems 
with two strange quarks (S) are more compact than the 
systems with only one strange quark (like E), and the lat- 
ter systems are more compact than the nucleon in charge 
density distributions. The statement becomes more clear 
when we look at the magnetic squared radii, where we 
have now 0.7 fm^ for N, 0.5-0.6 fm^ for E, and 0.3-0.4 
fm'^ for S. Note that, for {rf.j), the magnitudes hold also 
for the neutral particles. 

Our results for {r^) are compared in Table IVTl with 
the experimental results for p, n and E~. With the ex- 
ception for the neutron, to be discussed later, our results 
are consistent with the experimental values. We compare 
also our results with the estimates from Ref. [3l| . In that 
work the octet baryon electric charge radii were extrap- 
olated from the quenched lattice QCD data [63] to the 
physical point using chiral perturbation theory, includ- 
ing corrections for both the finite volume and quenched 
effects. Aside from A and E° which were not estimated 
(and n to be discussed latter), the results of Ref. [31| 
are in agreement with our results. We recall that A (like 
E°) has not been included in the calibration of our model, 
since the properties of A were not calculated in the lattice 
simulations in Ref. [52] , on which our parametrization is 
based. Note however, our result for the A electric charge 
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squared radius is also very small (0.08 fm^), although 
larger than the resuh of Ref. [65l|, 0.003-0.017 fm^. We 
recall again that the electric charge form factors of neu- 
tral particles are difficult to simulate in lattice QCD, due 
to the cancellation among the contributions from the dif- 
ferent flavor quarks, which should cancel out exactly at 
Q^ = 0. Our results can also be compared with other 
estimates presented in the literature [ll|, [12, [2J, I72H79J . 
Using the decompositions of Eqs. ([71| and ([75]) . we 
separate also the effects of the valence quarks and those 
of the pion cloud. The results decomposed for {r^^) are 



(rl), (r|)^ (r|) Ref [31] 



Exp. 



presented in the second and third columns in Table I VII 
Focusing first on the charged baryons, the valence quarks 
give dominant contributions, although the pion cloud ef- 
fects can be important and as large as 33% for the S+. 
Since the effective contributions of the valence quarks 
were also estimated in Ref. [31j, we compare {r^)^ di- 
rectly with those results (third column) in Table IVIIII 
Although the meson cloud contributions in Ref. [31] are 
not the same as ours, and also kaon cloud contributions 
are included (but amounted to small contributions), it 
can be illustrative to compare our results of {r%)^ with 
the estimates of Wang et al 1311 . Our results are surpris- 
ingly close to those of Ref. \3V\ , with a notable exception 
for S+. In this case our estimate for {r%) is 0.2 fm^, to 
be compared with that of Ref. JU, -0.061 ± 0.045 fm^, 
with opposite sign. In general, we conclude that our esti- 
mates for the valence quark contributions are consistent 
with those of Ref. ^. 

In Table rvm we present also the results for the bare and 
pion cloud contributions for the magnetic squared radii. 
Although the physical radii can differ appreciably from 
the lattice extracted radii without chiral extrapolations 
(see figures from CTTU)) . lattice results can give us an idea 
on the magnitude of the valence quark contributions. For 
this, we compare our results with those of the lattice in 
for the lowest pion mass, rriTr = 306 MeV, and 
for TO^ = 354 MeV. It is still interesting to 
notice that (f|/)j^ approach the lattice extracted values 
when the number of the valence strange quarks increases. 
The result for E is closer to that of the lattice than that 
of the nucleon, and the results for S is nearly equal to 
that of the lattice, when the standard deviation is taken 
into account. This is not unexpected, since lattice QCD 
calculations use nearly the physical strange quark mass 
value. These facts also give us some confidence on our 
model to estimate the valence quark contributions both 
in the physical and lattice regimes. 

We now comment on the neutron electric charge 
squared radius. Our result, —0.029 fm^, differs appre- 
ciably from the experimental value —0.12 fm^. This de- 
viation is a consequence of our global fit, and the absence 
of high accuracy data until recently for Geu in the low Q^ 
region [8^ . The result of our fit is presented in Fig. [TTl 
and magnified in Fig. [T51 for the region Q^ < 1.5 GeV^. 
As already explained, our fit includes the neutron physi- 
cal data, but also includes the octet baryon lattice data. 
In particular, the fit also includes the lattice data for the 




p 0.704 0.066 0.770 0.685(66) 0.769(8) [66] 

n -0.098 0.070 -0.029 -0.158(33) -0.1161(22) [51] 

A -0.0081 0.091 0.082 0.010(9) 

S+ 0.503 0.214 0.717 0.749(72) 

E° -0.0020 0.049 0.048 

E" 0.519 0.113 0.632 0.657(58) 0.61(15) [70] 

H" 0.090 -0.0036 0.087 0.082(29) 

H~ 0.404 0.019 0.423 0.502(47) 



TABLE VI: Electric charge squared radii of the octet 
baryons. All values are in fm^. The subindexes indicate re- 
spectively, bare (&) and pion cloud (vr) contributions, as de- 
fined in the text. The experimental value presented for the 
proton is an average of Refs. [6a - [6a ]. An estimate of the pro- 
ton electric charge squared radius using muonic hydrogen [691 ] , 
gives 8% less than the result presented in the table. 



neutron form factors Geu and Gmu- However, since we 
have reduced the impact of the lattice data for the neutral 
particles to compensate the inaccuracy of the data, the 
lattice constraints to the valence quark contributions are 
not very strong. As a consequence, the strongest con- 
straint for the neutron electric form factor comes from 
the very accurate, high Q^ physical neutron data. As one 
can observe in Fig. 1151 our model result deviates from the 
data in the region Q^ < 0.5 GeV^. In the figure the result 
associated with the valence quark contributions is closer 
to the data than the full result (dressed). We then con- 
clude that the deviation of our result from the physical 
data is a consequence of our estimate for the pion cloud 
contributions, which are poorly constrained by the low 
Q^ neutron physical data. In short, the poor description 
of Geu is a consequence of our ambition to describe si- 
multaneously the lattice and physical data. It would be 
improved if we concentrated only on the physical data 
in the low Q^ region, as done in the other studies |28j . 
The new generation of experiments for Geu with high 
accuracy [80] will be important to clarify the impact of 
the different effects in the neutron data, and can con- 
strain quark models such as the one used in the present 
work. Nevertheless, our analysis shows the importance 
and sensitiveness of the pion cloud effects on the neutron 
and the neutral baryons in general. The importance of 
the pion cloud effects is also manifested in the results of 
the electric form factors of S", A and S°. 



VII. CONCLUSIONS 

In the present work we have applied a covariant spec- 
tator quark model to study the valence quark structure 
of the octet baryons. Combining the contributions of the 
valence quarks (covariant model) with a phenomenolog- 
ical, covariant parametrization for the pion cloud con- 
tributions, we have described the electromagnetic form 
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(rji), (rji)^ (rli) Ref. [65] Ref. [52] 



p 


0.679 


0.059 


0.738 


0.470(48) 


0.40(8) 


n 


0.714 


0.0034 


0.718 


0.478(50) 


0.46(11) 


A 


0.544 


-0.242 


0.302 


0.347(24) 




E+ 


0.383 


0.146 


0.530 


0.466(42) 


0.36(8) 


E" 


0.365 


0.103 


0.468 


0.423(38) 




S~ 


0.407 


0.199 


0.606 


0.483(49) 


0.37(8) 


=•0 


0.370 


-0.0010 


0.368 


0.384(22) 


0.32(2) 


H~ 


0.295 


0.045 


0.340 


0.336(18) 


0.29(4) 



TABLE VII: Magnetic dipole squared radii of the octet 
baryons. All values are in fm^. The experimental results 
are (»"|/) — 0.733 ± 0.096 fm^ for the proton, and (tm) — 
0.767 ± 0.123 fm^ for the neutron [fH. More recent re- 
sults for the proton are (rfj) = 0.604 ± 0.026 fm^ [H, and 



(rli) = 0.752 ± 0.035 fm^ [69]. While the lattice resuHs from 
Ref. [63] correspond to the pion mass m^ — 306 MeV, those 
from Ref. [52l correspond to m-^ = 354 MeV. 





(-1), 


Ref [31] 


Ref. [65] 


Ref. [52] 


p 


0.704 


0.746(69) 


0.452(53) 


0.41(4) 


n 


-0.098 


-0.097(31) 


0.029(41) 




A 


-0.0081 


0.026(9) 


0.025(7) 




E+ 


0.503 


0.820(75) 


0.503(54) 


0.44(3) 


E° 


-0.0020 




0.046(15) 




E- 


0.519 


0.586(56) 


0.410(37) 


0.357(16) 


■=0 


0.090 


0.113(27) 


0.062(22) 




2" 


0.404 


0.471(46) 


0.389(23) 


0.326(6) 



TABLE VIII: Valence quark contributions for the octet 
baryon electric charge squared radii, compared with the va- 
lence quark contributions from Ref. [sj, and lattice QCD 
simulations [53 . I65( ] . While the lattice results from Ref. 65] 
correspond to the the pion mass ?n,r ~ 306 MeV, those from 
Ref. [52] correspond to m^ — 354 MeV. 



factors of the octet baryons in the lattice and physi- 
cal regimes. The octet baryon systems are described in 
a simplified quark model with the S-state configuration 
for the quark-diquark relative motion. The simplicity 
of the wave functions allows us to calibrate the momen- 
tum dependence of the octet baryon form factors using 
four independent range parameters. The model has been 
calibrated using lattice QCD data for the octet baryons 
(272 data points), the physical nucleon form factor data 
(202 data points), and the physical octet baryon mag- 
netic moment data (5 data points). Overall, we have 
used 6 parameters for the valence quark structure, and 7 
parameters for the pion cloud effects. We have derived a 
parametrization for the octet baryon wave functions de- 
pending on the baryon flavors. The parametrized wave 
functions are very useful for future studies of the electro- 
magnetic reactions. 

Our results suggest a dominance of the valence quark 



effects, but the inclusion of the pion cloud effects im- 
proves the global description of the octet baryon electro- 
magnetic form factor data. The pion cloud effects are 
of the order of 15% for the nucleon in the low Q^ re- 
gion, but decrease in the higher Q^ region. Surprisingly, 
the pion cloud effects are still important for the nucleon 
electric form factors in the region of 5-10 GeV^ with the 
contributions of 5-8%. This feature is a consequence of 
the value of A2 (Aj = 1.54 GeV^) for the pion cloud 
parametrization, which gives a slow falloff for the Pauli 
form factor F2. 

We have a good description of the nucleon data (x^ per 
data point —1-9), but for the description of the neutron it 
is slightly worse (x^ per data point of 1.9 and 2.3 for the 
electric and magnetic form factors, respectively). This is 
a con seq uence of the high accuracy of the recent Gmu 



data |6j], and also our ambition to describe simultane- 
ously the lattice and physical data . The large number of 
the lattice data points reduces the impact of the nucleon 
physical data. Also the inaccuracy of the neutron lattice 
data (especially differs from the neutron charge zero at 
Q- — 0) makes it difficult to constrain those parameters 
of the model directly related to the nucleon system. To 
improve the description of the nucleon system, it will be 
useful to perform a more systematic and detailed study 
for the nucleon form factors. This can be done using a 
more complete lattice QCD database including the re- 
sults from Refs. ISli - ISa ]. and also the results of the new 
generation of nucleon form factor data in the low Q^ re- 
gion [so, |83, [Sa] , which determine Ge and Gm simultane- 
ously instead of the ratio^ Ge/Gm- In this work we have 
restricted to use the lattice data from Ref. [53 , since one 
of our goals is the overall description of the octet baryon 
electromagnetic structure, and we prefer to avoid pos- 
sible inconsistencies arising from different lattice QCD 
simulations. 

Our results for the neutron form factors require some 
discussions. Our global fit has turned up to give a 
poor description for the neutron electric form factor for 
Q^ < 0.5 GeV^ (electric charge squared radius of —0.029 
fm^). Although Gsn is poorly constrained by the lattice 
data, it is complemented by the physical data. Our result 
for the neutron electric squared radius differs apprecia- 
bly from the experimental value, although our estimate of 
the bare valence quark core is closer to the experimental 
value, and also to that of the valence quark contribution 
from Ref. [3l[ . We interpret this as a consequence of the 
impact of the high Q^ region data (very precise in gen- 
eral), which reduce the impact of the low Q^ region data 
(less accurate) in the global fit. As already discussed, 
a more precise calibration of the model for the neutron 
will require a more detailed analysis of the nucleon sys- 



® That will require the combination of the polarization trans- 
fer measurements and cross section measurements (Rosenbluth 
method). 
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cQ^ 0.06 




o"o.04 



Q^ (GeV^) 



FIG. 15: Neutron electric form factor in the low Q^ region (mag- 
nified from Fig. [TTJ. The lines are for the total result (solid line) 
and the bare result (dashed line). 



tern, and this will be possible when more precise data 
as in Ref. [80j become available. More accurate lattice 
QCD data for the neutron can also be very useful. Nev- 
ertheless, our results show that the pion cloud effects are 
indeed important and influential for the neutron data. 

Overall, we have a very good global description of both 
the lattice QCD and physical data for the charged parti- 
cles p, S^ and E!~. For the analysis of neutral particles, 
one must be careful, because they have been less con- 
strained by the lattice data (n and E!°), or not at all con- 
strained (A and E"). But we can conclude that the pion 
cloud effects are very important for the electromagnetic 
form factors Geb of neutral particles. 

Concerning the size of the octet baryons, we conclude 
that, as expected, the systems with two strange quarks 
are more compact than those of the one strange quark, 
and the latter is more compact than the nucleon sys- 
tem with no strange quark. The estimates of the octet 
baryon sizes (electric charge squared and magnetic dipole 
squared radii) due to the valence quarks (bare core) are 
consistent with this conclusion, and are also in agreement 
with the independent estimates |31i] . The inclusion of the 
pion cloud effects in the electric charge squared radii, in 
general, makes the final results approach to the experi- 
mental results, and also to the values estimated based on 
the lattice results extrapolated using chiral perturbation 
theory ^. 

The agreement between our estimates with the octet 
baryon electromagnetic form factor data, and also those 
of the radii, gives us some confidence about the applica- 
bility of our model to estimate the valence quark effects 
both in the physical and lattice regimes. An application 
of the present model for the octet to decuplet electro- 
magnetic transitions is in progress by generalizing the 
study made for the jN — ^ A transition J18l - |20l | . For this 
purpose, we can use the octet baryon wave functions ob- 
tained in this work, and those of the decuplet obtained 
in Ref. H within the same formalism 18911. 
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Appendix A: Brief review of nucleon 
electromagnetic form factors 

Here we briefly review the results of Ref. [l| for the 
nucleon electromagnetic form factors. We consider the 
wave function from Eq. ^ , corresponding to the nucleon 
case (see Ref. [i[ for details). Below, M^ represents the 
nucleon mass. 

In the following we use 



^"^ = E=n(^)4:(A), 



(Al) 



where ep^ is the diquark polarization vectors of the final 
(P+) and initial (P-) states in the fixed-axis representa- 
tion [331 . The direct calculation gives [l,, il8i, ,35j : 



A"/3 = _ 



^a/9 



papP 
P+P- 



(A2) 



/ P P 

r, f P" + - P" 



Pi 



P+P- 



M\ 



,t2 



-Pt 



with 



M+M. 



P+ ■ P^{M+M_ + P+ ■ P_) ' 



(A3) 



where M_ and A/+ are the masses of the initial and final 
states. In the present case M+ = M_ = M^f. 



J 



Equation (|11|) leads to the transition current [1|, 
3 



N 



-B{Q') X u{P+) 



JiT 



32- 



i<j'"'qa 



2Mn 



(A4) 



;7a75 



J3T + J4 



2Mn 



757;3A"^|«(P_), 



where A"'^ is given by Eq. (jA2p . and 



BiQ' 



ijN{P+,k)^JjN{P-,k). 



(A5) 



The coefficients, ji and J2, can be written in the SU(2) 
sector, 



^/i+(Q') + ^/i-(Q')t3, 



(A6) 
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and 



3{i+2) = -^TjJ^Tj = -/l+ - g/l-T-3. (A7) 

Note that the expressions (jA6l) and (jA7[) . when ap- 

/ 1 ^ 



plied to the nucleon isospin state x 







X- 



1 



are equivalent to Eqs. ((T7)) and (fTS)) . 



With some spin algebra we obtain: 



Fi 


m 


= 2S(Q') 




[- 


1 1 




3 UP^ + Q2 


F2 


(Q^) 


= IsiQ') 


X 


{- 


1 1 


3 4M2, + Q2 



and 



(A8) 



[16M2,j3 + (4A4 - 3Q2)j4] 



To extend these results to the octet baryon B, we use 
Mb for its mass, and make the replacements: 



Ji 



J2 



J3 



J4 



Mat 
Mb^s 



(AlO) 
(All) 
(A12) 
(A13) 



The expressions associated with ji and j3 are determined 
directly from the respective definitions. For J2 and J4, we 
need to take into account that the quark current Eq. (fTi)) 
is written in terms of the nucleon mass, which leads to 



2Mjv Mn 2Mb 



(AM) 



■A,S 



That is why the coefficients J2 ' are modified by the 
(A9) factor ^j^. Finally, in Eq. (jA5P we replace -07v by V'b- 
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